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Abstract 

The classical and relativistic Hamilton-Jacobi approach is applied to the 
one-dimensional homogeneous potential, V(q) = aq n , where a and n are 
continuously varying parameters. In the non-relativistic case, the exact an- 
alytical solution is determined in terms of a, n and the total energy E. It 
is also shown that the non-linear equation of motion can be linearized by 
constructing a hyper geometric differential equation for the inverse problem 
t(q). A variable transformation reducing the general problem to that one of 
a particle subjected to a linear force is also established. For any value of n, it 
leads to a simple harmonic oscillator if E > 0, an "anti-oscillator" if E < 0, 
or a free particle if E = 0. However, such a reduction is not possible in the 
relativistic case. For a bounded relativistic motion, the first order correction 
to the period is determined for any value of n. For n » 1, it is found that the 
correction is just twice that one deduced for the simple harmonic oscillator 
(n = 2), and does not depend on the specific value of n. 

*rose@dfte. ufrn.br 
Tj anilo@dfte.ufrn.br 
^limaj as@astro.iag.usp.br 

1 



I. INTRODUCTION 



The Hamilton- Jacobi (HJ) equation is a powerful method either to the relativistic and 
classical framework, and also plays a prominent role in quantum mechanics as a special 
route for a comprehension of the Schrodinger equation [1,2]. The associated action variable 
method, which is a slight modification of the HJ approach, is also quite useful for determining 
frequencies and energies of periodic systems. Even for the general relativity theory, the 
importance of the HJ approach has long been recognized by many authors (see, for instance, 
[3] and Refs. therein). Indeed, one of the most elegant methods for describing geodesies and 
orbits of test particles in Schwarzchild and Kerr spacetimes, as well as for any stationary 
gravitational configuration, is provided by the relativistic HJ equation [3-5]. 

In the nonrelativistic domain, the standard applications of the HJ theory are the har- 
monic oscillator, the Kepler problem, and charged particles moving in electro-magnetic fields. 
Recently, the method has also been applied to the rocket problem [6]. Pars [7] studied the 
motion of a classical particle in a plane under central attraction derived from the potential, 
V(r) = ar n , for some particular values of n. However, to the best of our knowledge, if a 
and n are continuously varying parameters the general solution has not been obtained even 
for the one- dimensional case. 

In this letter we discuss an analytical solution for the power law potential both for the 
classical and relativistic case. As we shall see, the HJ equation for a particle moving under 
the action of a one-dimensional potential, V(q) = aq n , has a general and unified solution in 
terms of hypergeometric functions. The period of the motion as a function of the power index 
n and of the a parameter are easily determined from the general solution. For completeness, 
we also show that such a problem may also be exactly solved starting from the nonlinear 
equation of motion by employing a slightly modified Euler-Lagrange approach. By changing 
the coordinate and adopting an auxiliary time, we show that the motion of the particle in 
such a potential may be reduced to the problem of a particle under the action of a linear 
force: an oscillator if E > 0, an "anti-oscillator" if E < or a free-particle if E = 0, where 
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E is the total energy of the particle. Classically, these results hold regardless of the value 
of n. Such a reduction cannot be implemented in the relativistic case. However, the general 
problem may be reduced to an integral, which is a natural extension of the earlier Synge's 
treatment for the relativistic harmonic oscillator [8]. 



II. HAMILTON- J ACOBI APPROACH 

The classical Lagrangian for a particle subject to the one-dimensional potential, V(q) = 
aq n , reads 

X 

L{q,q,t) = -mq 2 - aq n , (1) 

where an overdot means total derivative. The parameters a, n, vary continuously, and m, q 
are, respectively, the mass of the particle and its generalized coordinate. The Hamiltonian 
is given by (jp is the canonical momentum) 

H(q,p) = ^- + aq n = E, (2) 

and since it does not depend explicitly on the time, it represents a conserved quantity which 
is the energy of the system. The Hamilton- Jacobi equation for the Hamiltonian (2) assume 
the following form [1] 

1 fdS\ 2 „ dS 

where S is the Hamilton's principal function. Since the explicit dependence of S on time is 
involved only in the last term, the variables can be separated. Following standard lines, the 
solution is assumed to have the form 

S(q,E,t)=W(q,E)-Et, (4) 

where E, the constant of integration, has been identified with the total energy. With this 
choice, the time is readily eliminated from (3), and Hamilton's principal function S becomes 
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S = j ^2mE - 2maq n dq - Et . (5) 

It should be noticed that if the energy E is positive, the constant a appearing in the 
potential may assume either negative or positive values, however, if E < 0, the allowed 
values of a are necessarily negative. Since H = E, such considerations also follows naturally 
from the positivity of the kinetic energy 

K = ]- mq 2 = E - aq n . (6) 

The second integration constant is a consequence of the Jacobi transformation equation 



b-—= r™~ r dq (7) 

' " dE \l2\E\Jo ^1 - 1 ' 



The root of \E\ arises naturally if one takes for negative values of E, E* = —E, and compute 
the partial derivative with respect to E. In order to integrate (7) we make the change of 
variable u = (a/E)q n . Thus we have that 

It is worth noticing that the integral appearing on the right hand side of the above expression 
is the incomplete Beta function [9] 



B u (a,b) = [ U y^il-yf^dy 
Jo 



for a = 1/n and b = 1/2. Without loss of generality we may take q(t = 0) = 0, so that (3 = 
in (8). When the motion is finite, the particle moves back and forth (oscillatory motion), 
and the period T can be easily determined. We see from (6) that at the turning points 
(q = 0) the amplitude of the oscillation is A = (E/a) 1 l n ) and, therefore, u max = 1, is the 
maximum value in the upper limit of the integral (8). The period of the motion is given by 
T = At m ax and, taking these considerations into account, we find that equation (8) yields 



T = 2J^-A 1 - 1 / n B(l/n,l/2), (9) 



an 2 
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where B(l/n, 1/2) = y^T (j^j /T ^ + |) is the complete Beta function. The period of this 
oscillator is clearly dependent on the amplitude A, unless n = 2. As should be expected, for 
n = 2 equation (9) reduces to 



T = 2tT, 



m 
2a ' 



which is the well known result for the simple harmonic oscillator. 

It should be recalled that the incomplete Beta function is related to the hypergeometric 

Gaussian F(a, b; c; u) by the following identity: B u (a, b) = a~ 1 u a F(a, 1 — b; 1 + a; u). It thus 
follows that the result (8) may be expressed as 

f +^&-(^ + 1 ^")' (10) 

As a check, we notice that for n = 2 the above expression can be written as 

u(t + P)=zF(±l2*) , (11) 



where we have defined the variable z = \Jct/ E q and u = 2n/T is the frequency of the 
oscillatory motion. Since F(l/2, 1/2; 3/2; z 2 ) = ^arcsin^: [9]; by returning to the old 
variable q, the above equation can be recast as 



/ 2E 

9(*) = \/ o sinw(t + /3), (12) 



moo 2 



which describes the motion of a simple harmonic oscillator [1]. 

III. MODIFIED EULER-LAGRANGE APPROACH 

We present here another method of solution for the problem outlined in the introduction. 
In principle, by considering that the HJ approach leads to the hypergeometric function for 
t(q) as given by (10), it should be possible to obtain the differential hypergeometric equation 
starting directly from the equation of motion. The importance of this method is two-fold: 
the initial conditions remains arbitrary and the linearization procedure is a guarantee that 



5 



we get the complete solution for the nonlinear differential equation governing the motion of 
the particle. 

The equation of motion (Euler-Lagrange equation) for the Lagrangian given by (1) is 

mq + naq n ~ l = . (13) 

Instead to solve the equation of motion above, we will consider an equivalent equation. By 
multiplying it by q, and inserting the definition of the total energy of the particle given by 
(6), we obtain: 

n .o nE ^ 
2 m 

It is worth mentioning that, if some few identifications are made, (14) is the same differential 
equation describing the evolution of the scale factor in the standard Friedman-Robertson- 
Walker (FRW) cosmological model [10]. 
A first integral of (14) is 

t = ™(i-!Lf). (15) 



m V E 

It should be noted that the signals for the total energy E and for the constant a must be 
properly chosen. So, if E > the constant a can have any signal, but if E < we must 
have a negative a. The integration of (15) is achieved with the introduction of an auxiliary 
variable defined by u = (a/E)q n . With the aid of this transformation the first integral is 
written as 

^ = nS« 1 - 1 /»(l- u )V2 > (16) 

(-lib 



Fm /E\ 1/2 - 1/n 



where 



m \a J 

At this stage we could integrate (16), as it was done in Sec. (II). Instead, on may consider 
the inverse problem, i.e., the solutions for a differential equation expressing the time as a 

6 



function of the auxiliary variable u(q). This inversion is readily done through (16), which 
provides the first derivative of t(u): 

dt 1 



— = —u 1/n -\l-u)- 1/2 , (17) 
du nB v ; K J 

whereas the second derivative may be put in the form 

«1_ )<* +[C _( 8+1M *_0. (18) 

where c = 1 — 1/ n and a — l/2 — l/n. This is Gauss's hypergeometric differential equation. 
Its general solution consists of a linear combination of two linearly independent solutions 
[12]: 

t(u) = CiF(a, 0; c; u) + C 2 v}~ c F(l + a - c, 1 - c; 2 - c; u) , (19) 

where F is the hypergeometric function, and C±, C 2 are the arbitrary constants of integration 
to be determined by the boundary conditions. It is a property of the hypergeometric function 
that if any of the first two parameters is zero, then the series terminates (i.e., F(a, 0; c; u) = 
1). It thus follows that the solution (19) reduces to 

' = c 1 + ^f) 1/ V(ii;i + i;|<f), (20) 

where we have substituted the value of the parameters a and c. The case n = 2 gives 

q(t) = J- sinw(t-fo), (21) 
V a 

where the constant uj — I/C2 and to = Ci remain to be determined. 



IV. PARAMETRIC SOLUTIONS AND PERIODICITY 

As we have seen, if E ^ 0, the general solution t(g) cannot be inverted to obtain explicitly 
q(t). In such cases, parametric solutions are usually more enlightening. Some simplicity is 
achieved when we replace the pair (t,q) — > (t,Q), where the "conformal time" r and the 
new coordinate Q are defined by 



dt = q{r)dT , (22) 

and 

Q = qS . (23) 

Under the above transformations, the expressions for the Lagrangian and energy E, given 
by (1) and (2), become 

l = ~tQ~ 2 Q' 2 — a Q~ 2 , (24) 

E = 2 -^Q-^ + aQ- 2 , (25) 

where a prime denotes derivative with respect to r. As one may check, the Euler-Lagrange 
equations are now given by 

QQ"-Q' 2 -f^ = o , (26) 

and inserting Q' 2 from the energy conservation law, we obtain 

O--g£) O = 0. (27) 

This is an interesting result. The above equation describes the motion of a classical particle 
under the action of a linear force. If E < (E > 0) the force is of restoring (repulsive) type, 
while the motion of a free particle corresponds to E — 0. The general solution of (27) is 

Q = % sin yfe(iUT + 5) , (28) 

where e = —2E/m, uj = n/2 is the frequency of the oscillator (anti-oscillator) and Q , 5 
are integration constants. Since the energy E may be negative, the "phase" 8 may assume 
complex values. The constant Qo can be determined using the energy equation written in 



terms of the conformal time. It turns out that Qo = ^ — 2a/m. Substituting this into (28), 
and returning to the original coordinate q, we have 

?W = (- m /2a)V»(Ei^l)" 2/ ", (29) 
where we have put 5 = in (28). 



V. RELATIVISTIC HAMILTON- J ACOBI APPROACH 



Let us now consider the relativistic treatment for the homogeneous potential. For a 
single-particle the relativistic Lagrangean is 



L = -mc 2 Jl-^-aq n . (30) 



The Hamiltonian function for the particle is given by the general formula 

,dL 



H = q—-L = Jp 2 c 2 + m 2 c 4 + aq n , (31) 
oq v 

where p is the canonical momentum. The Hamilton- Jacobi equation for S is obtained by 
replacing in (31) p by dS/dq and H by —(dS/dt): 

C2 (f) 2 "(f + <*f + mV = 0. (32) 
As before, a solution for (32) can be found separating the variables in the form 

S(q,E,t) = W(q,E)-Et, (33) 

where the integration constant E is again the total energy. With the above choice one finds 
from (32) 

S=\J y/(E - aq n ) 2 - m 2 cHq - Et , (34) 

which should be compared with its classical version (5). The second integration constant 
arises out from the transformation equation ((3 = |^) and we have finally 

t + P = ±[ , E - aqU dg. (35) 

cJ ^{E -aq n ) 2 -m 2 c A 

Without loss of generality we take q(t = 0) = 0, so that (3 = in (35). We next suppose that 
the potential is a symmetric function about the origin, that is, we consider only even values 
for n. Then the motion will be limited between V(—A) and V(A) where the amplitude is 
now given by A = [(E — mc 2 )/^ 1 /", and, from (35), the period T is determined by 



T 



4 r' 
c Jo 



dq 



(36) 



'i /_mc 2 _V 

\E-aq n J 

Writing the total energy as E = mc 2 (l + e), and considering that the potential energy is 
small compared to the rest mass energy mc 2 , relativistic corrections to second order in the 
period T are 

dq 



nn 4 rA 
T — - 

c 



> =(l + ^4"-<f)) 

o ^2n(A n - q n ) 4 



(37) 



where k = a/mc 2 . Introducing the variable change y = (q/A) n the above integral can be 
rewritten as 

4 



T 



O A ,l + 2 

f l / . _/ In , 3/t/l 2 /■! ,2-1, i 

/ y 1 » j (l -y) { *>dy + =■ / y ( n >(1 -y)2dy 

Jo c«v2k Jo 



(38) 



cAt ~ 1 uV2k Jo " cu^/2k Jo 

whose values are tabulated in terms of the complete Beta function [9]. Substituting for k 
we write the period as 

3 aA n 



T = 2\ 



2m 



an* 



A x - nt2 B{\ln, 1/2) 



1 + 



2n 



(39) 



8 mc 2 \n + 2 

Observe that for n = 2 (relativistic harmonic oscillator), the above equation reduces to 

T = 2^\l + 3 -^\ , (40) 
V 2a I 8mc 2 \ v ; 

where the term (3aA 2 )/8mc 2 is the first order relativistic correction earlier obtained by 

Synge [8] following a different approach (see Appendix). 

The general form (39) is an interesting one. The term multiplying the square bracket is 

just the nonrelativistic period T for the bounded homogeneous potential as given by (9). 

The remaining term represents the relativistic correction 

A^WJrM (41) 
T 8mc 2 \n + 2J ' v ; 

and, although it depends of n, we see that for n ^> 1 it saturates around (3a A n ) / Amc 2 which 

is twice the relativistic correction for the harmonic oscillator. 

In conclusion, we have shown that the classical version of homogeneous potential problem 

can completely be solved by the HJ approach. In the relativistic case we not succeed in 
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obtaining an analytic solution, but the earlier Lagrangian Synge's treatment for the harmonic 
oscillator is readily generalized using the HJ transformation. Finally, we stress that although 
important from their own right, such problems may also have interest for describing some 
excited states appearing in quantum chromodynamics which have been usually described by 
the relativistic oscillator [13]. 
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VI. APPENDIX: REMARKS ON SYNGE'S INTEGRAL 

The relativistic harmonic oscillator (n = 2) was long ago discussed by Synge [8]. In this 
appendix, we show how the exact Synge's integral form for the period of the motion can be 
generalized for a power-law index. 

For an arbitrary n, the period is given by (see (36)) 

T. = ±r , ^ „ ■ (42, 



c Jo 



' 1 _ ( >"c 2 \ " 
\E-aq n ) 

Now, by writing E = mc 2 + aA n one finds 



T = 2 /— [ A [i + 2x 2 (i-(gM)")]^ (43) 

V aA n Jo ^1 - ( q /A) n + X 2 (l - {q/AYf ' 

where x 2 — aA n /(2mc 2 ). 

Introducing the variable <p by q = A sin 2 /™ ip, the above integral becomes 



Tn = A\[^A l - nl2 ^ Sin2/n_1 * (1 + C ° s2 V?) ^ 



i-n/2 H sm ' <f[L-^Zx cos yjdy? 
an 2 Jo \/l + x 2 cos 2 (p 



and, for n = 2 we have 



t 2 = r (1 i 2x2c ; s2 ^ , (45) 

V 2a Jo v 1 + x cos 2 cp 



which is the Synge's integral with a slight different notation (our a = k 2 /2, and A = a in 
Synge's work [8]). Finally, by expanding (44) in powers of x-, the relativistic correction for 
the period as given by (39) is recovered. For n=2 it reduces to Synge's expression. 
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